H — T phase diagram of the bi-dimensional Ising model with exchange and dipolar 

interactions 
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We explore the equilibrium properties of a two-dimensional Ising spin model with short-range ex- 
change and long-range dipolar interactions as a function of the applied magnetic field H. The model 
is studied through extensive Monte Carlo simulations that show the existence of many modulated 
phases with long range orientational order for a wide range of fields. These phases are characterized 
by different wave vectors that change discontinuously with the magnetic field. We provide numerical 
evidence supporting the existence of first order transitions between these phases. At higher fields 
our results suggest a Kosterliz-Thouless scenario for the transition from a bubble to a ferromagnetic 
phase. 



I. INTRODUCTION 

Thin magnetic films have been the subject of intense 
attention over the last two decades [iHll- Most studies 
have been motivated mainly by the technological appli- 
cations of these structures But, they also faced sta- 
tistical physicists with the challenge of trying to answer 
many foundational questions regarding the role of micro- 
scopic interactions in the macroscopic behavior of a large 
system. 

These quasi two-dimensional structures show a large 
variety of ordering effects including formation of striped 
states, reorientation transitions, bubbles formation in 
presence of magnetic fields and hysteresis At the 

origins of these phenomena is the competition between a 
short-ranged interaction favoring local order and a long- 
range interaction frustrating it on larger spatial scales. 
The role of the long-range interaction is to avoid the 
global phase separation favored by the short-ranged in- 
teraction promoting, instead, a state of phase separation 
at mesoscopic or nano-scales. Then, it is not, in gen- 
eral, a small perturbation, [8l| but must be considered as 
precisely as possible. 

From a computational point of view, this means that 
the frustrating interaction has to be accounted for by 
involving all the lattice sites in the computation, which 
in turn limits the actual system size that can be handled 
in Monte Carlo simulations. On the other hand, to obtain 
exact results on multi-scale, multi-interaction systems is 
extremely difficult, so that simulations are often the only 
source of information. 

Model Hamiltonians taking into account short-ranged 
exchange ferromagnetic and long-range dipolar anti- 
ferromagnetic interactions have been used to reproduce 
many of the elemental features observed in experiments 
of magnetic systems @. Unfortunately, and despite the 
obvious relevance from the experimental point of view 



of the presence of an external magnetic field, most of 
the numerical studies so far have concentrated their at- 
tention on the zero magnetic field case {H = 0). This 
is in part because of the already very rich and complex 
phenomenology obtained by tuning the strengths of the 
exchange and the dipolar interactions, but also because 
of the almost prohibitive computational cost of the sim- 
ulations, even for moderated lattice sizes. 

To fill this gap, we use extensive Monte Carlo sim- 
ulations to determine the role of an external magnetic 
field in the thermodynamical properties of quasi two- 
dimensional magnetic systems. We present results for 
systems where exchange and dipolar interactions are 
comparable and where the anisotropy contribution to the 
Hamiltonian is very large. 

The work is organized as follows: In section [ll] we 
present the model and review some of its properties. In 
section IIIII we give details about the Monte Carlo sim- 
ulations and discuss the motivation for the parameters 
used and its connections with previous reports in the lit- 
erature. Then, in section HVl we present and discuss our 
results. This section is organized in three parts, we first 
present and analyze the H — T phase diagram of the 
model, then we provide some insight on the ground state 
structure of the different phases, and finally we charac- 
terize the transitions between these phases. Finally, in 
section |V] the conclusions of the work appear. 



II. THE MODEL 

We consider a square lattice of Ising spins oriented per- 
pendicularly to the plane of the lattice and interacting 
through the dimensionless Hamiltonian 
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where Si = ±1 is the value of the spin at site i. The first 
sum runs over all pairs of nearest neighbor spins and the 
second over all pair of spins in the lattice. The discrete- 
ness of Si is consistent with infinite or very large magnetic 
anisotropy The parameter S = Je/Jd stands for the 
ratio between the strength of the exchange and dipolar 
interactions, Je and Jd respectively. H is the magnetic 
field intensity (in units of Jd) and rij is the distance, 
measured in crystal units, between sites i and j. 

This model, but in zero external magnetic field, has 
been extensively studied^lO 14]. For example, it is now 
well understood that in a wide range of values of S its 
ground state consists in stripes of anti-parallel spins with 
a width that increases with S Once the tempera- 

ture is turned on, the situation becomes more complex 
and, in a (5 — T phase diagram, one can recognize a zool- 
ogy of phases, stripes of different widths, paramagnetic 
phases, tetragonal, smectic, nematic, and others fu. IT3j. 
Roughly speaking, at zero field the system presents a first 
order phase transition between a low temperature phase 
of stripes and a high temperature tetragonal phase with 
broken translational and rotational symmetry. It was 
also shown [ll| that for a narrow window around (5 = 4 
the model develops a nematic phase where the system has 
short range positional order but long range orientational 
order. 

On the other hand, in ref. [15| Garel and Doniach 
study analytically the H — T phase diagram of a contin- 
uous Landau-like model with dipolar interactions. They 
conclude that the H — T plane is characterized by three 
different phases: stripes, bubbles and ferromagnetic. 
Their analysis also suggests a scenario with Fluctuation 
Induced First Order Transitions (FIF0T)^6,] between 
the phases, or a second-order melting of the Kostcrlitz- 
Thouless (KT)[l3| type for the bubble- ferromagnetic 
transition. While some of this phenomenology is con- 
firmed by our simulations, we will show below that the 
phase diagram resulting from the Hamiltonian ([1]) is even 
richer. 

Numerical simulations using Langevin Dynamics [3~ 
[20j on similar Landau-like models seem to support the 
general picture described in [l^. In particular, in ref. 
|19| the author studies the behavior of the system under 
external magnetic field, but focus his attention mainly on 
the role of metastable configurations, the presence of hys- 
teresis loops and memory effects. Therefore, the predic- 
tions of ref. [15,] are still waiting for conclusive numerical 
support. 

For the particular case of the Hamiltonian ([IJ , the cor- 
rectness of the predictions of Garel and Doniach [Ih] is 
even less clear. While at first one expects that the cor- 
respondence between the standard ferromagnetic Ising 
model and the continuous 0^ model persists even in the 
presence of the dipolar term, the existence of commensu- 
ration effects, typical of striped patterns in discrete Ising 
systems, may alter this intuition. For example, the au- 
thors of reference [TJ] studied the H — T phase diagram 
of Hamiltonian ^ using Montecarlo simulations and 



found no evidence for the transition to a bubble phase, 
suggested a continuous character for a stripe-tetragonal 
boundary and reported some unexpected jumps in the 
magnetization versus temperature curves. 

In this sense our work revisits these previous simu- 
lations looking with more attention to the effect of the 
magnetic field at low temperatures. Some of the results 
already seen in jl4] are confirmed and, we think, analyzed 
in more detail and from a differeiit perspective. Some re- 
sults support the predictions of [la], and others, to our 
knowledge, are new, and enrich the already complex phe- 
nomenology of these systems. 



III. SIMULATION 

We centered our analysis on the value (5 = 4 which 
corresponds to a zero field ground state of perfect al- 
ternating stripes of width h = 2. So, for the smaller 
system sizes considered we have 8 periods of modulated 
stripes. In all cases the size of the system L was properly 
commensurate with the period of the H = modulated 
phase. Arlett et al. [l3| used values of 6 between 6 and 
8, having stripes of width h = A and 6 respectively. So, 
for the system sizes they consider 4 or at most 6 periods 
of the modulated structures are present. As we will dis- 
cuss below this makes difficult to interpret some of the 
consequences of the presence of H . 

This value of 5 = 4 is representative for proved first 
order stripes-tetragonal transition in if = but it is 
also known to be of the order of real mag netic- frustrated 
systems seen in experimental works [10']. Some connec- 
tions between experimental systems and values of rela- 
tive strengths of interactions in theoretical models can 
be found in reference [2lj. More recently, Carubelli et 
al. qualitatively reproduced detailed measurements 
of magnetic changes of samples of Fe/Ni/Cu(001) [l^ 
by means of a Heisenberg-spins model, very similar to 
Hamiltonian ([ij, using a value of ^ = 6. 

To build the phase diagram, the system is first initial- 
ized in the equilibrium configuration at a fixed tempera- 
ture and zero magnetic field. To guarantee equilibration 
the magnetic field is increased very slowly 10^* < AiJ < 
10~^ and for each [H, T) point, we let the system relax 
for ti = 10^ Montecarlo steps (mcs) using a Metropo- 
lis dynamics. Once equilibrated, the system evolves over 
other t2 = 10^ mcs to measure the physical quantities 
of interest. We impose periodic boundary conditions to 
limit finite size effects and explore different values of tem- 
perature and field for linear system sizes up to L = 48. 
To account for the long-range interactions we implement 
the Ewald Summation Technique [23>] adapted to the par- 
ticular case of the magnetic dipolar potential j2J] . 
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IV. RESULTS AND DISCUSSION 

A. Phase Diagram 

The main result of this work is shown in figure [TJ This 
is the H — T phase diagram of the model represented by 
the Hamiltonian ([1]). 
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FIG. 1: Phases diagram for a system of L=32 considering the 
anharmonic zone; void circles and squares transition lines are 
first order. The critical temperature for the transition to the 
tetragonal phase at 77 = is shown with a colored square, the 
high temperature zone is represented with a slight shadow. 

Four different zones are well defined in the diagram. 
For low values of temperature and external magnetic field 
the system is in a oriented modulated phase of perfect 
stripes characterized by a wave vector k = (0,7r/2) and 
zero magnetization. Increasing the magnetic field, new 
modulated phases, characterized by new wave vectors, 
and non-zero magnetization appear. These new phases, 
keep the orientational order but are characterized by sev- 
eral wave vectors (therefore we call them anharmonic 
phases) that depend on the magnetic field. The prop- 
erties of these phases and the location of the transitions 
suffer from strong finite size and commensuration effects, 
so, in the diagram we represented only one zone that, for 
the system size considered, contains all the anharmonic 
structures. Similar phases were already predicted within 
a mean field scenario for an Ising model with competing 
interactions Jq and Ji between nearest and next nearest 
neighbours in one direction of a cubic lattice (ANNNI 
model) m. 

For still larger values of H we find a phase without 
orientational order (bubble). Finally, increasing further 
the magnetic field the system becomes completely mag- 
netized (ferromagnetic phase). At low H, and close to 
the stripe to tetragonal transition the combination be- 
tween thermal fluctuations, commensuration and finite 
size effects, and the excitations due to the magnetic field 
makes the analysis of the phase diagram too difficult. So, 
in this zone, the structure of the phase diagram is still 



unknown, and we shadow this zone in figure [T] to caution 
the reader about this. 

Now, to fix the ideas, let us concentrate our attention 
on the results for one temperature. We define, following 
p^ . the so called n/2 rotational symmetry-breaking (SB) 
parameter: 

^7= : 1 (2) 

+ Uh 

where Uv [nh) is the number of vertical (horizontal) 
bonds between nearest neighbors anti-aligned spins. This 
parameter takes the value I in a perfectly ordered stripe 
state while it equals zero for any phase with tt/2 rota- 
tional symmetry. 

In figure[2]we represent the evolution of as a function 
of H for T = 1.2 in a system with iV = 32 x 32 spins. The 
zooms show typical configurations for the corresponding 
values of rj. As can be seen, abrupt jumps separate clear 
plateaus of rj at three different values of the magnetic 
field, H ~ 0.84, H ~ 1.34 and H ~ 2.40. Each plateau 
reflects an underlying symmetry of the system. 




FIG. 2: SB parameter in a system of L = 32 as a function of 
the field for T = 1.2, the spots are some typical configurations. 

A deeper understanding of the phase diagram and spe- 
cially on the character of the jumps separating the dif- 
ferent plateaus is obtained analyzing figure [3] where the 
magnetization, the magnetic susceptibility and the sus- 
ceptibility associated to the rotational order parameter 
77 {2Txri =< 7f > — < rj >^) are plotted. Increas- 
ing from zero the external magnetic field, the rotational 
symmetry-breaking parameter rj and the magnetization 
show various plateaus separated by abrupt jumps (see fig. 
[3^). These jumps result from the discrete properties of 
the lattice where the model is defined. Discrete changes 
in the field are required to change from one stable struc- 
ture of stripes to another. 
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The existence of these jumps is also clearly reflected 
in both susceptibilities (see figures [5}3 and [5J;). Three 
different peaks are well defined in the magnetic and the 
orientational susceptibility at the same transition points 
where the orientational order parameter and the magne- 
tization jump. 

For H > 3, the magnetization starts to growth lin- 
early with H but the rotational symmetry-breaking pa- 
rameter is zero. The system is in the so called bubble 
phase already predicted by Garel and Doniach [15] for 
the Ginzburg-Landau model with dipolar interaction. Fi- 
nally at very high fields {H > 5) the system is completely 
magnetized, m = 1 and 77 = 0. 
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FIG. 3: Evolution under increasing magnetic field for a system 
of L—32: a) SB parameter and magnetization for T = 1.2, 
b) magnetic susceptibility for T — 1.2 and c) SB parameter 
associated susceptibiUty for T = 1.2. 



The jumps in the order parameter and the peaks in the 
susceptibilities suggest the existence of different thermo- 
dynamic phases at each plateau of rj. To characterize 
the properties of these phases we look at the form of the 
structure factor, S{k), in each plateau : 
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FIG. 4: Evolution of structure factor in a system of 1/ = 32 
for increasing magnetic field. 

Figure H] shows the structure factor of the system for 
different values of H . Each plot is obtained by the av- 
erage of 5000 equilibrium configurations. At very low 
magnetic field, the system is characterized by a peak at 
one wave vector k — (0,7r/2). Increasing H new peaks 
appear in S{k). First, with component k = (0,fcy ^ 0), 
still signaling the presence of orientational order in one 
direction. This change in the form of the structure factor 
is not evident a priory. One may, for instance, expect that 
the external magnetic field unbalances the number of up- 
down spins creating defects that breaks the orientational 
order. Our results suggest a different scenario, where 
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if properly equilibrated at low temperatures, new struc- 
tures, without evident defects, keep the orientational long 
range order of the original ground state structures. Then, 
at higher magnetic fields, (see in the figure H = 2.52) 
S{k) becomes symmetric in both axis, the system looses 
the orientational order and reaches the bubble phase. Fi- 
nally the magnetization saturates and only the peak at 
k = (0, 0), survives. 

Figure [S] shows the contribution of the three principal 
wave vectors k = (0,fc*) characterizing the evolution of 
the system configurations with the magnetic field. Ini- 
tially, the perfect stripes phase is characterized, as we 
know, by a wave vector k = (0,7r/2). At « 0.84 a 
new wave vector k = (0,77r/16) dominates the system, 
still indicating the presence of oriented stripes. Increas- 
ing further the magnetic field, at iJ « 1.34, S{k) changes 
again, and k = (0, Stt/S). The sudden rise and decay of 
each wave vectors reflects again the abrupt changes in 
the symmetry of the system. 




FIG. 5: Main harmonic contributions to the equilibrium con- 
figuration structure factor for increasing H in a system of 
L = 32. Other ky contributions remains always under 0.2. 

We also calculated the directed spatial correlation 
functions for the system 
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which reveal interesting information about the equilib- 
rium states. In particular, we tried to fit the numerical 
data with a function of the form 

C{r)=Ae~'^cos{kr-^) + Br-°'+D (4) 




FIG. 6: Fits of correlation numerical data with the function 
(Q; correlation best fits for two particular values of apphed 
field 



that has been proposed for the approximated continuum 
model[26, 27]. Figure |6] shows the corresponding fits for 
averaged equilibrium configurations at two values of H. 

From these fits we can gain information about the de- 
pendence with H of the correlation length of the modu- 
lated domains (^), the main wave vector of the phase (fc) 
and the power law strength (a) respectively. In particu- 
lar, we can see in figure [7] the behavior of fc as a function 
of H . The plateaus in k coincide with the principal wave 
vectors (see figure [S]) characterizing the different stripe 
structures. 



o 




FIG. 7: 

L = 32. 



Best fitted k value as a function of // in a system of 



To our knowledge these new anharmonic phases have 
not been predicted before in a model with dipolar inter- 
actions. They are absent in the continuous model, where 
the effect of the magnetic field in the striped phase, is 
considered assuming that below the bubble phase the 
stripes persist in an increasing magnetized background 
[T5| . They are present in the ANNNI model, but differ- 
ently from [U the ANNNI model is anisotropic by con- 



6 



struction. 

On the other hand, in the phase diagram resulting from 
the simulations in |14l | the orientational order parameter 
changes continuously from a finite value to zero at a given 
field (see figure 7 in that reference) . The reasons for these 
differences in the phase diagrams are not clear. We are 
tempted to think that looking at the dependence of 77 for 



lower values of the temperature the authors in ref. 14 1 
could find similar jumps and phases. Of course, having a 
large S and hence larger stripe widths the anharmonicity 
properties of their structures may be hidden by strong 
finite size effects. 



B. Ground state analysis 

To study what kind of structures are responsible of the 
anharmonic phases, we tested the energy of a large num- 
ber of configurations of alternating S = —1 and S — 1 
stripes. The width of the S = —I stripes was varied from 
1 to 2 while the width of 5' = 1 stripes was varied from 
to L, as it is expected for the striped configurations in the 
presence of a field H > 0. Thus, borrowing the notation 
from ref. \2§\ we denoted as one configuration with 
stripes of width 2 against the field and stripes of width 
4 in the field direction, repeated periodically. 

The energies of these configurations are represented in 
figure E] as a function oi H. At H = 0, the ground state 
of the system corresponds to the h2 phase. By increasing 
H the system reaches a critical field Ha, where perfect 
stripes becomes energetically unfavorable with respect to 
the anharmonic configuration (ft,2322, ky — Tir/W). For 
larger fields, a new anharmonic configuration becomes 
the ground state (/i32, ky — Svr/S). Further increasing 
H the situation repeats with the appearance of new an- 
harmonic states. How many of these anharmonic con- 
figurations may appear depends strongly on temperature 
and commensuration effects. The corresponding ground 
state energies of the system, considering only these an- 
harmonic configurations is represented in figure [S] with a 
continuous line. This line corresponds to the lower energy 
curve obtained from the superposition of the energies of 
the different configurations as a function of H. 

One may wonder whether these are finite size effects, 
and a non-orientatied ground state structure may dom- 
inate the behavior of the infinite system al low H. To 
test our predictions, this analysis was repeated for dif- 
ferent system sizes, N ^ L x L. Figure IH] suggests that 
independently of the system size, the first critical field 
Ha always appear in the low field region where the per- 
fect stripes become unstable. This value defines a zone 
in which anharmonic structures establish, mainly in the 
form of ft,2223 or h23 configurations depending on com- 
mensuration effects. 

The transition between anharmonic and bubbles 
phases remains around H — 2A for system sizes up to 
L = 48, this have been used to draw a schematic broken 
line in figure [S] Since the bubble phase establishes be- 
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FIG. 8: Energy as function of in a system of L = 32. 
Different anharmonic configurations at T = and the ground 
state energy (continuous line) are shown. Ha is the critical 
field at which perfect stripes are lost. 
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FIG. 9: Critical fields changing configurations at T = ver- 
sus system size. Ha and ferro fines are bolded for clarity but 
vertical lines were the only region explored and no interpola- 
tion is obvious. An schematic bubble line is also drawed (see 
the text). 



couse entropic effects, this is lakely to be valid for large 
system sizes. In all tested cases the critical field remains 
well below this schematic transition, supporting the ex- 
istence of the anharmonic phases obtained for L — 32 
in the thermodynamic limit, and giving rise to a rather 
wide anharmonic zone. 



C. Phase Transitions 

Unfortunately the computational cost associated with 
the presence of long range interactions and the commen- 
suration effects in this kind of systems, prevent us from 
doing a proper finite size scaling analysis to define the 
character of the transitions. Instead, we focused our at- 
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tention in systems of sizes L = 32 and L = 40 and study 
the histograms of the energy and the order parameter. 



Evidence for First Order Phase Transition 

The jumps in the susceptibihties and the discontinu- 
ities in 7] in figure [3] aheady suggest the first order charac- 
ter of the transitions between the different orientational 
phases, and from the last anharmonic phase to the bub- 
ble phase. However, a stronger evidence is given in fig- 
ure [ini These histograms were calculated for systems of 

= 40 X 40 spins, sampling lO^mcs after relaxation for 
each value of H and considering lO'^ values of energy and 
V- 
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FIG. 10: Histograms of the orientational order parame- 
ter and energy for fields around the transitions involving 
striped phases: a) and b) for harmonic-anharmonic, c) and 
d) anharmonic-anharmonic, e) and f) anharmonic-bubbles. 



For the three transitions considered, the figure shows 
that, increasing 7?, the histograms of the order param- 
eter and energy evolve from unimodal functions at low 
magnetic fields, to a two peak shape structure, that dis- 
appears at higher magnetic fields giving rise to the new 
thermodynamic phase. For the particular case of the 
anharmonic-anharmonic transition (c-d), the difference 
in energies between the two structures is so small that the 
histograms for the energy appear always as unimodal. 

On the other hand, one must note that while in the 
first two transitions, the peak in P{E) moves from high 
to low energies, in the anharmonic to bubble transition 
it moves from low to high energies. In this transition, 
the system looses the orientational order and therefore E 
increases. This is compensate by the presence of strong 
entropic effects that, in this more disordered structure, 
dominate the equilibrium state of the system. 

It is relevant for the definition of the anharmonic to 
bubble transition the appearance at high fields of a non- 
zero correlation length ^ for the modulated domains (see 
figure [TT|) . Fitting the spatial correlations in the bubble 
phase with expression Q we obtain the expected inverse 
proportionality of ^ with the applied magnetic field [2^ . 
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FIG. 11: Best fit of correlations by expression Correlation 
length of modulated domains as a function of H in a system 
of L = 32, for H below the bubbles region ^ > L. 



Evidence for a KosterUz-Thouless transition 

Figure [12] shows some views of the domain structure 
of the system close to the bubble-ferromagnetic transi- 
tion. They suggest that increasing H the bubble phase 
dilutes in a ferromagnetic environment. This support 
the predictions in J^] where the authors proved that 
within a Ginzburg-Landau approximation, dislocation of 
the bubbles structure may lead to a second-order melting 
transition of the Kosterliz-Thouless type. The continu- 
ous change of energy and magnetization (see figure [3^) 
and the saturation of the response functions close to this 
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transition (see in figure fT3l zooms of the magnetic suscep- 
tibility and the specific heat close to this transition) also 
support these predictions. 



H=6.0 






FIG. 12: Some field values involved in the transition from 
bubble to ferromagnetic phases. Typical configurations for 
an L=48 system. 
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FIG. 13: Field values involved in the transition from bubble 
to ferromagnetic phases. Specific heat and magnetic suscep- 
tibility versus field in a system of L = 32. 

One last indication in favor of this scenario, comes from 
the spatial correlation functions of the system. In figure 
[14] we show the value of a obtained by the fits of the cor- 
relation functions with expression The sudden rise 
of a close to if 4.5 is also consistent with a Kosterliz- 
Thouless transition. 



CONCLUSIONS 



We present evidence supporting the idea that all, but 
the bubble to ferromagnetic, are first order transitions. 
This is also in agreement with analytical results that pre- 
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FIG. 14: Best fit of correlations by expression Q. Power 
exponent as a function of Ji" in a system of L = 32. 



dieted that the stripes to bubbles transition is of the Bra- 
zovskii [bi] type. On the other hand, close to the bubbles 
to ferro transition, our simulations show the existence of 
a continuous order parameter, the saturation of the re- 
sponse functions and algebraically decaying spatial cor- 
relations, supporting all, a Kosterliz-Thouless scenario. 

Finally, it is worth to note the interesting parallelism 
between these anharmonic phases and the hybrid states 
found for Hamiltonian ([1]) at zero field in reference [29| . 
There, through Mean Field calculations, the authors 
suggested a possible interpretation of nematic phases 
as a competition between striped structures of differ- 
ent widths. Moreover, they found Kosterliz-Thouless 
features in the transition between striped and nematic 
phases. To clarify these issues, and to completely de- 
fine the phase diagram ([T|) more accurate simulations are 
expected close to the H — Q critical temperature. 



We developed extensive numerical simulations to char- 
acterized the phase diagram of the model given by ([IJ. 
This Hamiltonian presents at very low field and temper- 
ature a phase of symmetric stripes and zero magnetiza- 
tion. Increasing the field, new thermodynamical phases 
appear, still with orientational order but with non-zero 
magnetization and characterized by different wave vec- 
tors. As far as we know, the existence of these ther- 
modynamic phases have not being proposed before for 
systems with dipolar interactions. For larger values of 
H, the system enters into the bubble phase loosing the 
orientational order. Then, at larger fields, the system 
becomes fully magnetized. 
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